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Section 1 


Introduction 

1*1 Purpose 

The DSAS Illumination and Moon Conflict programs vill provide the times 
during an orbit when the DSAS (Digital Solar Aspect Sensor) vill record 
the direct rays of the sun, and the periods of time when the horizon 
scanners vill come in conflict vith the moon. The DSAS Illumination 
Program makes use of an orbit tape (or epoch time and orbital elements) 
in addition to an ephemeris tape containing positions of the sun and 
moon. The Moon Conflict Program makes use of the same ephemeris tape 
vith sun and moon positions, but uses only epoch time and orbital 
elements for the satellite positions. These programs vere designed for 
the TIROS or ITOS series spacecraft but may be utilized by any spacecraft 
vith similar sensors. 
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1.2 General Characteristics 


The DBAS Illumination Program uses orbital information, either from an 
0RB1 tape or Brouwer mean orbital elements, and ephemeris positions of 
the sun. The Moon Conflict Program uses Brower mean orbital elements 
with a Brower orbit generator and ephemeris positions of the moon. These 
programs are not entirely restricted to the TIROS (ITOS) series of satellites. 

With proper adjustment of mounting angles and testing times, these programs 
could be used on satellites with configurations and sensor equipment similar 
to TIROS (ITOS). The only restrictions in these programs are that the 
satellite spin axis pitch angle must be near zero degrees in the DSAS Illumination 
Program and that the spin axis have a nominal attitude of in the 

Moon Conflict Program. Further information on the configuration of the 
TIROS (ITOS) spacecraft and mission objectives can be found in the attitude 
determination system documentation generated for support of these satellites. 

This can be obtained from the GSFC Computer Program Libraxy, filed under 
Program Humber 000198 (TADS). The DSAS Illumination Program ("SUHILL") 
is filed under G00392; the Moon Conflict Program ( "MORCON" ) under G00393. 

2. DSAS Illumination Program 
2.1 Analysis 

ftie test for illumination of the Digital Solar Aspect Sensor (DSAS) will 
be made wing the known sun position and the satellite orbital position 
at a predetermined time. This test will be performed at any increment of 
time provided the time step is in whole minutes. Under normal operating 
conditions the satellite attitude is assumed nominal, that is, Yav (<{>) = 

Roll ($) ■ Pitch (0) ** 0°. Allowance is made for a variation in the nominal 
attitude by inputting the current values for $, ♦ and 0. 

The satellite position in inertial coordinates is calculated from the pre- 
dicted orbital position using the equations: 
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A —A —A 

n = R x V 

IrxvT 

t = n x r 

where R- the inertial satellite position vector 
V= the inertial satellite velocity vector 
r= unit position vector in inertial coordinates 
n= unit vector normal to R and V in inertial coordinates 
i* unit vector normal to N and R in inertial coordinates 

If the satellite attitude is ^=$=0=0°, then the satellite spin axis attitude 

in inertial coordinates is calculated using: 

A A 

s=n 

* 

where s = spin axis vector in inertial coordinates (S_, S_, S.). If the satellite 
attitude is not nominal, then the right ascension («) and declination (6) of the 
spin axis must be input. In this case, the spin axis attitude in Inertial coordinates 
is obtained from: 

f S x = cos * cos 6 
® = j 8y»8in « cos 6 
l S z * sin 6 

To determine if the DBAS will be illuminated, the sun's vector must be known in 

the satellite coordinate system. Given the position vector of the satellite (from 

orbital data) and the sun's vector to the earth (from an ephemeris tape) , both 

at the desired point in time, the sun vector to the satellite is obtained by 

the equation: 

P* Q - R~ 

|Q - R| 
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where P » unit satellite to sun vector 
q * earth to sun vector 

(see Figure l) 

If it is assumed that the satellite does not have a nominal attitude, that 
is, 4 + 0° and/or 4 4 0° and/or 6 + 0°, then the spin axis vector § resolved 
in the rotating orbit frame by: 

e A 

S r « r« s 
S t - t . S 

A A 

Sg* n . s 

where (S r , S t , S n ) are the orbit frame coordinates of the satellite spin axis. 

The satellite pitch angle 0 is automatically maintained near 0° by the satellite. 

4 and 4 must be solved for in terms of S r , S^, Sq. This can be accomplished by 
two consecutive right-handed rotations of the tody frame vector (0, 0, 1) through 
the angle (-4) and then (-4)* 4 and 4 are solved using: 

with quadrant check 

4*sin -1 (s r ) 

A 

Having determined the angles 4» 4 and 0, the sun vector P is then resolved in 
the body frame by first resolving it in the orbit frame: 

ft-i'KVVrV. 

V” ■ p ~ n x p z* a y p r n z p z 

Then , transform P Into the hod, frame using three right-handed rotations through ♦ , 
4 and 6: 

Pi’ - P r 

» tt f • 

P 2 «Pfc cos 4 + P q sin f P 1 B " p 3 sin ♦ + P 1 cos 4 

l n l . . 

P3 * -P t sin 4 ♦ P n cos 4 (l) ?2 * ^2 '2; 

II • • 

P3 • Pj cos 4 ♦ Pj^ sin 4 
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SATELLITE IN SUNLIGHT TEST 


SATELLITE 



RE = Earth's Radius 

R = Radius Vector of Satellite Orbital Position 

-A 

Q = Earth to Sun Vector 

P = Unit Satellite to Sun Vector 

Satellite will be in Sunlight if ft - > 0 


FIGURE - 1 



(3) 


Pj=>P^" cos 0 + P2 sin 0 
„ •» 

t 2 “ sin 0 ♦ Pp cos 0 
P 3 » P 3 " 

where (P]_, P2, P3) are the coordinates of the sun vector in the body frame. 

(see Figure 2) 

If the satellite has a nominal attitude, that is, ** $ = 0“ 0°, then equations 
1, 2 and 3 will reduce to: 

p l sP r 

P 2“ P t 

P 3* P n 

To determine if the DSAS will record the sun, several tests must be made. The 

first is a test to see if the satellite is in sunlight .Using the satellite 

position vector R and the position coordinate of the sun vector in the orbit 

frame (P r ) , one can calculate the angles y + 6 using* 

s'* arc sin ) 

|R I 


where R« • earth's radius in kilometers 
$* arc cos (-P r ) 

(see Figure 1) 

The satellite will be in sunlight if B~y>0. 

The next tests involve the DSAS body frame sensing plane vectors D and E. D 
is the unit normal vector and E is the unit centerline vector. These vectors are 
resolved in the body frame by: 

Dj* cos B 
D2 m sin B 

d 3 « 0 

Ej» sin B 
Eg® cos B 

e 3 » 0 
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BODY ORIENTATION ANGLES 


7,V (OUTWARD RADIAL) 



EULER SEQUENCE * YAW ABOUT r to F=r, 2 ' , T f i 0° <T < 360° 

* ROLL 0 ABOUT V to I s , 2V = 2? , Ti -90° < 0 < 90° 

• PITCH 9 ABOUT T to 1, 2, 3=3"? -90° < 9 < 90° 

NOTE: For Nominal Mission Mode Attitude, 

T, "2, T coincide with"?, "t, n, andY= 0=9 = 0 

FIGURE - 2 
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vhere B» the mounting angle of the DSAS measured from year axis. (D 3 and are 0 

since the sensor is mounted 90° to the pitch (spin) axis.) (See Figure 3.) A 
test must he made to determine if the sun's rays vill strike the sensing side 
of the DSAS. The angle between the DSAS sensing centerline E and the sun 
vector P is found by: 

AEP * arc cos (E*P) 

A A 

where AEP * the angle between E and P. If AEP<90°, then the sun's rays will strike 
the sensing side of the DSAS. 

The next test will determine if the sun's rays fall within the field of view 

A 

of the DSAS. The angle between the spin axis and the sun vector P is found 
by: 

ASAP* arc cos (P 3 ) 

A 

where ASAP* angle between spin axis and P. If the lower viewing angle of the 
DSAS is greater than ASAP and ASAP is less than the upper viewing angle* then 
the sun's rays are within the field of view of the DSAS as measured from the 
positive spin axis. 

When all the preceeding tests are passed, the only remaining test is to determine 
exactly when the sun is recorded. In this test the angle between the sun vector 

A A 

P and the DSAS sensing plane normal D must be calculated from: 

D0T*J).P 

where cosine of the angle between D and P . This test will not be passed 
until the sun ' 8 rays cross perpendicular to D and the cosine changes in sign. 

'I'M? hen is the interval of time when the DSAS will record the sun pulse. 

2.2 Setup 

2.2.1 Resources 

1. Card reader 

2. Three tape drives (two if orbital elements are used) 

3. Three tapes — system tape, sun and moon ephemeris tape, 

and orbit tape (omit orbit tape if orbital elements are used.) 
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DSAS MOUNTING AND VIEWING ANGLES 


+ 2 
ROLL 



+ 3 



For TIROS-M Viewing Readout Range for DSAS is 128® 

FIGURE - 3 
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k. Card, deck - JCL and data cards. 

2.2.2 JCL before data (see Figure 4) 

// JOB card (use 2 minutes CPU and I/O time estimate) 

//A EXEC LOADER, PARM»'EP=MAIN' ,REGI0N=250K 

/ /SYSLIN DD UNIT-(2U00-9,,DEFER),LABEL»(l,BLP),DISP*(0LD,KEEP), 

//DCB=(RECFM=FB,LRECL»80 .BLKSIZE-3200 ,DEN-2 ) ,VOL=SER»33 1 »62J 

This JCL is for the 9-track system tape containing both the DSAS 
Illumination Program and the Moon Conflict Program. The DSAS Illumination 
Program is located in File 1 of the system tape, the Moon Conflict Program 
in File 2. As of January 1973 there exist two duplicate system tapes. The 
tape numbers, either of which is to be inserted after "SEP*", are 33462 J 
and 33568J. 

//GO.FT28F001 DD DSN AME-GDRTAP,UNIT=2U00-9, LABEL- ( ,BLP) , 

/ / VOLUME-SER- ,DISP=(NEW,PASS)^dcB»( ,LR£CL-8276 , 

// BLKSIZE-82 80, DEN-2, RECFM-V) 

(FT28F001 describes the 9-track sun and moon ephemeris tape . The 
current (January 1973) tape number is 30278K, which is to be entered 
after "SER=".) 

/ /GO . GENTAP DD UNIT-2 400-7 , DISP- ( OLD .KEEP ) , LABEL- ( ,BLP ) , 

// DCB» ( , DEN-1 .BUFL-U212 .BLKSIZE-4212 ) ,VOLUME=SER» 

(GENTAP defines the 7-track orbital tape. Hie number is entered after 
"SER-" ). (if orbital elements are used, these two cards are to be 
removed) . 

//G0.FT09F001 DD SYSOUT»A,DCB=(RECFM»VBA,LRECL»137,BLKSIZE»7265) 
//GO.FTlOPOOl DD DUMMY 

(The FT09F001 and FT10F001 are the unit descriptions for the output 
units. The FT09 description prints the deta left-adjusted on the 
paper, the FT10 prints the data center-adjusted on the paper. Only 
one statement should be used, the other would then be set up aa 
DUMMY). 


//G0.DATA5 DD* 



2.2.3 Data Cards 


Card #1 

COL 

FORMAT 

VARIABLE 

DESCRIPTION 

1-2 

12 

N0RB1 

00:0 orbit tapes not used. Epoch 

time and orbital elements used. 

01 = orbit tape used. Epoch time 
and orbited elements flagged 
with 9's (see orbital elements 
and epoch time . ) 

Card #2 

COL 

FORMAT 

VARIABLE 

DESCRIPTION 

1-10 

F10.3 

HRDIFF 

See below 

11-15 

15 

ISMALL 

Time increment for testing illuminations 
( Normally 1 minute) 

16-20 

15 

IDELT 

Time jump in minutes after a conflict 
has occurred, this will minimize CPU 
in testing for the next illumination. 

2l»-30 

AU.A3 

ID6AT1 , IDSAT2 

7 alphanumeric characters designating 
the satellite name (s). 

HARDIFF»The difference between the input epoch tine and orbital elements and 
the desired epoch time . This difference is in hours. This variable 
will also update the elements to correspond to the desired epoch time. 
This variable is normally 0.000 with the decimal in column 7. By setting 
the variable to 0 . 000 , the time and elements will not be updated. This 
method should only be used to update the epoch time and orbital elements 
when the desired values are not available. 

Card #3 

COL 

FORMAT 

VARIABLE 

DESCRIPTION 

1-1* 

IU 

IEDATE(l) 

Epoch Year in It digits, e.g. 1972 

5-6 

12 

IEDATE(2) 

Epoch Month in 2 digits, e.g. 01 

7-8 

9 

12 

Leave blank 

IEDATE(3) 

Epoch Day in 2 digits, e.g. 15 

10-11 

12 

IEDATE(U) 

Epoch Hour in 2 digits, e.g. 12 

12-13 

12 

IEDATE(5) 

Epoch Minutes in 2 digits, e.g. 36 

lU-15 

12 

IEDATE(6) 

Epoch Seconds in 2 digits, e.g. 26 
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(Be careful to use ninutes and seconds. If orbit tape is used for input, 
be sure to make each digit a "9" e.g. IEDATE(l) =9999.) 


Card ik (Orbital Element s ) 


COL 

FORMAT 

VARIABLE 

DESCRIPTION 

1-10 

F10.2 

E(l) 

A, semi«maJor axis in kilometers 

11-20 

F10.5 

E(2) 

e, eccentricity (no units) 

21-30 

F10.3 

E(3) 

1, inclination in degrees 

31-U0 

F10.3 

E(U) 

mi mean anomaly in degrees 

fcl-50 

F10.3 

E(5) 

w, argument of perigee in degrees 

51-60 

F10.3 

E(6) 

6, fight ascension of ascending node in 
degrees 

(If orbit tape is 
9999999999). 

used as input. 

make each digit a "9", e.g. E(l)« 

Card #5 




COL 

FORMAT 

VARIABLE 

DESCRIPTION 

1-10 

F10.1 

B 

See below 

11-20 

F10.1 

ANGLOW 

M 

21-30 

F10.1 

ANGUPP 

tl 


B= Sun sensor mounting angle in degrees measured from Yaw axis to DSAS 

center line in direction of Roll axis, (reference spacecraft documentation). 

ANGLOW= Lower viewing angle limit of the sun sensor measured from the 
positive spin axis (equal to 26.5 degrees for ITOS spacecraft). 

ANGUPP* Upper viewing angle limit of the sun sensor measured from the 

positive spin axis (equal to 153.5 degrees for ITOS spacecraft). 


Card #6 

(Start 

date and time for program) 

COL 

FORMAT 

VARIABLE 

DESCRIPTION 

1-U 

IU 

ISTRT(l) 

Year 

5-8 

IH 

ISTRT(2) 

Month 

9-12 

IU 

ISTRT(3) 

Day 

13-16 

Ik 

ISTRT(U) 

Hours 

17-20 

Ik 

ISTRT(5) 

Minutes 

21-2 U 

Ik 

ISTRT(6) 

Seconds (The seconds are not used) 
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Card #7 (End date and tine for program) 


COL 

FORMAT 

VARIABLE 

DESCRIPTION 

1-4 

14 

IEHD(l) 

Year 

5-8 

14 

IEHD(2) 

Month 

9-12 

14 

IEHD(3) 

Day 

13-16 

14 

IEND(U) 

Hours 

17-20 

14 

IEHD(5) 

Minutes 

21-24 

14 

iehd(6) 

Seconds (the seconds are not used) 

Card #8 



COL 

FORMAT 

VARIABLE 

DESCRIPTION 

1-10 

F10.3 

PSI 

Yaw angle of spin axis in degrees 
(nominally 0.000) 

11-20 

F10.3 

PHI 

Roll angle of spin axis in degrees 
(nominally 0.000) 

21-30 

F10.3 

THETA 

Pitch angle of spin axis in degrees 
(nominally 0.000) 

31-40 

F10.3 

ALPHA 

Right ascension of spin axis in degrees 
(nominally 0.000) 

1*1-50 

F10.3 

DELTA 

Declination of the spin axis in degrees 
(nominally 0.000) 

/•(end 

of file) 




//(null) (Hot used if following JCL is used) 

2.2.4 JCL After Data 

If more than one copy of the listing is desired the Utility Program 
PATRICK may be used. The setup to invoke PATRICK is as follows: 

//STEP EXEC PGM=PATRICK,PARM*'9TN,001,001' ,REGION=100K 

// IH1 DD DSN= # . A. GO. FT09F001 ,DISP= ( OLD ,PASS ) 

(Uiis card references STEP A in the EXEC card at the beginning of 
JCL. If output is to be written to unit 10 instead of unit 09, then 
FT09 in the above step must be changed to FT10.) 
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NOTE: The FT09 definition most be changed if it is desired to write 

to the disk rather than SYSOUT=A. The cards should be: 

//GO.FT09F001 DD DSN*&&DUMl,DCM»(RECFMaVBA,LRECL»13T,BLKSIZE=T265), 

/ / DISP- ( NEW .PASS ) , SPACE- (TRK , ( 50 , 10 ) ) .UNIT-DISK 

(If FT10 is used, than change the statement to FT10F001) 

//0UT1 DD SYS0UT-A,SPACE-(CYL,(5,1)), 

// DCB»(RECFM«VBA,LPECL«137,BLKSIZE»7265 

These six cards may be repeated as many times as output is needed. Each 
set is an exact duplicate of the other. 

/• (End of file) 

II (Null) 

This completes the setup for the cards to run the DSAS Illumination 
program. 

2.3 Flowchart 

A flowchart of the general functions of the DSAS Illumination Program 
is given in Figure 5. This diagram shows the transformations and tests 
used to determine when DSAS (Sun Sensor) illumination has occurred. 



JCL DECK FOR DSAS ILLUMINATION PROGRAM 


//ZNSRMSUN JOB ( GH65 0160 IB ,P , G00392 , 002002 ) , AAA , MSGLEVEL = 1 
//A EXEC LOADER, PARM=' EP=MA IN' ,REGI0N=250K 

//SYSLIN DD UNIT = ( 2400-9 ,, DEFER), LA B EL= ( 1 , BL P ) , DI SP= < OL D, K EEP ) , 

// DCB=(RECFM=FB,LRECL=80 ,BLKS I Z E=3200, DEN=2 ) , V0L=SER=33462J 
/ /GO »FT28F001 DD D SN AME= GD RT AP , UN I T = 2400-9 , L A BE L= ( , BLP ) , 

// V0LUME=SER=30278K,DISP=< NEW, PASS) ,DCB=( ,LRECL=8276» 

// B LKSI ZE=8 28 0, DEN=2 , R ECFM = V ) 

/ /GO. GENT A P DD UN I T =2400- 7, DI S P= (OLD, K EEP ) , LABEL = ( , BLP ) , 

// PCB= ( ,DEN=I,BUFL=4212,BLKSIZE=42I2), VOLUME = S ER= 1234M 

//GO. F T09F00 1 DD SYSOUT=A ,DCB= ( RECFM=VBA , LR ECL = 1 37, BLKSI Z E=7 265 ) 

/ /GO .FT10F001 DD DUMMY 
//G0.DATA5 DD * 

01 


0.000 

1 50 

ITOS- D 



99999999 999999 





9999.99 

9.99999 

99 9. 999 

9.999 

999.994 

30.0 

26.5 


153 .5 



1972 10 15 

17 28 

00 




1972 10 18 

00 00 

00 




0.000 

0.000 


0. 000 

0.000 

0.000 

/* 






// 







FIGURE - 4 
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FLOW DIAGRAM OF DSAS ILLUMINATION PROGRAM 



/!rbit\ 


SUN & 
MOON 
EPHEM 


'mounting 

ANGLES, 

J START /END 
OF RUN 



I 


INITIAL 
ORB. DATA 
& SUN 
POSITION 



INITIALIZE 

PARAMETERS 

AND 

CONSTANTS 



CALC. SAT. 
POSITION IN 
INERTIAL 
COORDINATES 




SPIN AXIS POSITION 
IN INERTIAL COORDIN- 
ATES USING ORBITAL 
DATA 


RESOLVE 



SPIN AXIS POSITION 
IN INERTIAL COORDIN- 
ATES USING /, 6 


SPIN AXIS IN ROTATING ORBIT FRAME, 
CALCULATE SUN VECTOR TO SATELLITE 


CALC. SUN a 
VECTOR (P) 
IN ROTATING 
ORBIT FRAME 




FIGURE - 5A 
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FLOW DIAGRAM OF DSAS ILLUMINATION- PROGRAM(CONT'D) 



$ INTO BODY 
COORDINATES 


FRAME 

USING 
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2 . U Subroutines 


The following subroutines are used within the DSAS Illumination 


Program: 


BRRV 

XKEP 

BBRWR 

DJUL 

ALLOT 

MAD 

ALLOTZ 

RQUM2T 

BRWRl 

RDJPL 

BRWR2 

DAYS 

ELRV 

DPRB1 

ATARQ 

D2PRB1 

DMAD 

DELTAD 


The subroutines BRRV (main program) , BBRWR, ALLOT, ALLOTZ, BRWRl, 
BHWR2, ELRV, ATAHQ, DMAD and XKEP are used in the Brouwer mean orbital 
generator. The function of these subroutines is to provide orbital 
position parameters for the input satellite using epoch and orbital 
elements. Documentation for the Brouwer mean orbital generator can be 
found in "A Mutual Visibility Computer Program for Communication 
Satellites," (X-5UT-65-222) , by G. D. Repass and R. G. Chaplick, 

(May 1965). 

The subroutine DATS computes the number of days since the beginning 
of the input year for use in the orbit subroutine D2PRB1. 

The subroutines DJUL and MAD are used to calculate the Julian date 
for use in the Brouwer orbital generator and in the subroutines RDJPL 
and RQUM2T. The MAD subroutine performs the same function as the MOD 
function on the IBM 360 system. Subroutine DJUL is documented in the 
same document as the Brouwer orbital generator. 
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The subroutine RDJPL reads the JPL ephemerls tape for the inertial 
sun and moon position vectors at the input time. RQUM2T then rotates 
the sun and moon vectors to the true equator and equinox of that date. 

The subroutine DPRB1 reads the header record of a double precision 
seven track 0RB1 tape. D2PRB1 reads the double precision seven track 
0RB1 tape to obtain inertial radius and velocity vectors of the satellite 
position. 

Subroutine DELTAD was written for the DSAS Illumination and Moon 
Conflict Programs for the purpose of updating the time-keeping parameters 
within the main programs. The time variables used in DELTAD are: 

IDATE (1) ■ year IDATE (U) » hour 

IDATE (2) = month IDATE (5) = minute 

IDATE (3) * day IDATE (6) * second 

MINUTE » number of minutes to update time (5 digits maximum). Time vill 
be updated by any positive number of whole minutes up to five 
digits maximum. Updating for seconds was not needed and therefore 
not programmed. A general flow diagram is given in Figure 6. 

2.5 Output Description (see figure 7) 

A printed listing (multiple copies possible) is the primary output 
from the DSAS Illumination Program. If output on tape is desired, the 
JCL must be changed to meet this end. The output will consist of the 
following: 

1. The first page is a title: "Sun Sensor Illumination Test for * 
Satellite". 

*Ihe satellite name will be entered in this space from card #2. 
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FLOW DIAGRAM OF SUBROUTINE DELTA D 


' READ IN 
[DATE AND 
MINUTE , 


CONVERT 
MINUTE TO 
HOURS AND 
MINUTES 


UPDATE 

IDATE(4) 

AND 


'TDATE(5) 

'w 

yesT 


ID ATE ( 5) = 
IDATE(5)-60 
IDATE(4) = 
IDATE(4)+1 


UPDATE 


llDATE(3) AND IDATE(4) 
CORRECT FOR PROPER DAY 


UPDATE 

IDATE(2) 

AND 

IDATE(l) 


LEAP YEAR 


UPDATE 

IDATE(2) 

AND 

IDATE(3) 





2 


. lbe second page contains pertinent information: Sun Sensor mounting 
angle. Sensor viewing angle from positive spin axis, epoch time and 
orbital elements (9's are inserted if orbit tape is used), and start 
ing and ending date in year, month, day, hour, minute and second. 

3. The following pages usually contain tvo times during which an 
illumination occurred. These two times will be within one minute 
of each other. The format of the times are year, month, day, hour 
and minute. No seconds are written out. It is possible to obtain 
only one time, that is, when an illumination occurred on the minute. 
The output for this case will be one time and the message "Direct 
Illumination" . 
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Figure 7A 


22 



SUN SENSOR MOUNTING ANGLE 
(DEGREES) 

( 150 



Figure 7B 
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Section 3 


Moon Conflict Program 


3.1 Analysis 

The test for a moon conflict vith the earth horizon scanners will be 
made using the known moon position and the satellite orbital position 
at a predetermined time. Each test will be made at the time step 
increment beginning with the initial start time of the ran. The 
satellite attitude is assumed to be constant vith the Ysv (#)»Roll 
(0)»Pitch ( 0)«O°. These values cannot be changed unless the program 
itself is changed. 


The satellite spin axis vector is considered to be parallel vith orbit 
normal and is obtained from the radius and velocity vectors of the 


orbital position, R and V respectively. Hence, the unit spin axis is 
defined by: 

si-lil. 

IS * V| 

The unit vector from the satellite to the moon (M)can fce found using: 
A _ UTM - R 
|U™-F| 

where UTM* vector from earth to moon obtained from the ephemerls tape. 


The angle between the vector from the satellite to the moon and the 
vector from the satellite tangent to the moon can be determined. This 
angle will be called SIGMA(o) (see Figure 8) and will be calculated 
using: 

SIGMA “tan m 

(|M| 
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vhere R^radius of the moon in kilometers 

M*non-unit vector from the satellite to the moon 
The angle LAMBDA (X) between unit M and the positive spin axis unit 
vector can be defined by: 

LAMBDAscos** 1 (SA*M) 

(see Figure 8) 

Given the mounting angles of Scanner #1 and Scanner #2 as DELTA1 and 
DELTA2 respectively, a moon conflict will occur when the vector from 
the satellite tangent to the moon falls within the field of view of 
a scanner. This criterion will be met when: 

(1) DELTA1 - OMEGA - SIGMA * LAMBDA 4 DELTA1 + (MEGA + SIGMA 

or 

(2) DELTA2 - (MEGA - SIGMA 4 LAMBDA 4 DELTA2 + (MEGA + SIGMA 
where OMEGA=half viewing angle of the scanners 

Case 1 is a conflict test with the scanner with DELTA1 mounting 
angle while Case 2 is a conflict test with the scanner with DELTA2 
mounting angle. This latter case is illustrated in Figure 8 • noon 
petition 3. There will be a short period of time (a few hours) when 
the moon will lie between the viewing fields of both scanners but 
not in conflict with either one. This case is illustrated in Figure 
8 , moon position 2. 

Since the moon's orbital period is about twenty-eight days, it 
can be expected to cause conflicts every fourteen days. However, 
during one period of conflict, the moon will either be dark or 
only partially illuminated. To allow for this problem, the right 
ascensions of both the sun and moon are computed using: 
RAM00N=tan- 1 / uTM-(2) 

\utmTT); 
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MOON CONFLICT DIAGRAM 
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RASUN-tan" 3 / ) 
\UTS(l)/ 


where RAMOOH=right ascension of the moon 

RASUN*right ascension of the sun 

UTM(1)»X component of UTM (earth to moon vector obtained from 
the ephemeris tape) 

UTO(2)*Y component of UTM (earth to moon vector obtained from 
the ephemeris tape) 

UTS(l)*X component of the earth to sun vector obtained from 
the ephemeris tape. 

UTS(2)=Y component of the earth to sun vector obtained from the 
ephemeris tape. 

The difference betveen these right ascensions are computed from: 
DIFFERENCE® | RAMOON-RASUN | 

TOis difference is then converted to a MOD 180°. When the value obtained 
is greater than 90°, a true conflict will occur. If the difference is 
less than 90°, then the moon is only partially illuminated, and in 
the case of the TIROS-M (IT0S-1) satellite this was not considered 
to be a true conflict. Therefore , true conflict periods only occur 
once every twenty-eight days. 
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3.2 Setup 


3.2.1 Resources 

1. Card reader 

2. Two tape drives 

3. Two tapes (system tape and ephemeris tape) 

U. Card deck (JCL and data cards) 

3.2.2 JCL Before Data (see Figure 9) 

//JOB card (Use U minutes CPU and I/O time estimate) 

//B EXEC LOADER , PARM* ' EP=MAIN ’ ,REOION*250K 

(This is step name B for the Moon Conflict Program) 

//SYSLIN DD UNIT* (2^00-9,, DEFER), LABEL* ( 2, BLP ) ,DISP*( OLD, KEEP) , 

// DCB=(RECFM*FB,LRECL*80,BLKSIZE*3200,DEN*2) ,VOL*SER*33U62J 

(These cards reference the second file on the 9-track systems 
tape for the DSAS Illumination and Moon Conflict programs. 

The tape number to be inserted after SER-is either 33*»62J 
or 33568J (Note: these are duplicate tapes; either may be 
used) . 

//GO.FT28F001 DD DSNAME*GDRTAP ,UNIT»2U00-9 .LABEL* ( ,BLP) , 

// VOLUME=SER= , DISP= ( NEW , PASS ) ,DCB* ( LRECL-8276 , 

// BLKSIZE*8280,DEN*2,RECFM*V) 

(This is the 9-track ephemeris tape for the positions of 
the sun and moon. Currently the tape used is 30278K and is 
inserted after SER=.) 

//G0.FT09F001 DD SYSOUT-A ,DCB* ( RECFM-VBA , LRECL-137 ,BLKSIZE*726 5 ) 
/ /GO . FT10F001 DD DUMMY 
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(The JT09F001 and FT10F001 are the unit descriptions for the output 
units. The FT09 description prints the data left-adjusted on the 
paper, the FflO prints the data center-adjusted on the paper. Only- 
one statement should he used , the other cm he set up as DDHff • ) 
//G0.DATA5 DD* 

3.2.3 Data Cards 
Card #1 


COL 

FORMAT 

VARIABLE 

descriptioh 

1-10 

F10.3 

DELTA1 

Mounting angle of Scanner #1 from the 
positive spin axis measured in degrees 
(approximately 87.000 degrees). 

11-20 

F10.3 

DELTA2 

Mounting angle of Scanner #2 from the 
positive spin axis measured in degrees 
(approximately 93.000 degrees). 

21-30 

F10.3 

OMEGA 

Scanner half viewing angle measured in 
degrees (approximately 1.500 degrees). 

3h**h0 

Ah, A3 

ID6AT1 ,IDSAT2 

7 alphanumeric characters designating 
the satellite name(s). 

Card #2 




COL 

FORMAT 

VARIABLE 

DESCRIPTION 

1-U 

IU 

IEDATE(l) 

Epoch Year (e.g. 1972) 

5-6 

12 

IEDATE(2) 

Epoch Month (e.g. 12) 

7-8 

12 

IEDATE(3) 

Epoch Day (e.g. 12) 

9 

blank 



10-11 

12 

IEDATE(h) 

Epoch Hour (e.g. 01) 

12-13 

12 

IEDATE(5) 

Epoch Minutes (e.g. 25) 

lk-15 

12 

iedate(6) 

Epoch Seconds (e.g. 20) 
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Card #3 (Orbital Elements) 


COL 

FORMAT 

VARIABLE 

DESCRIPTION 

1-10 

F10.2 

E(l) 

a, semi-major axis in kilometers 

11-20 

F10.5 

E(2) 

e, eccentricity (no units) 

21-30 

F10.3 

E(3) 

i, inclination in degrees 

31-U0 

F10.3 

E(U) 

m, mean anomaly in degrees 

Ul-50 

F10.3 

E(5) 

v, argument of perigee in degrees 

51-60 

F10.3 

E(6) 

right ascension of the 
ascending mode in degrees 

Card #U 




COL 

FORMAT 

VARIABLE 

DESCRIPTION 

1-10 

F10.3 

HRDIFF 

See below 

11-15 

15 

MINUTE 

II II 

HRDIFF»The difference 

between 

the input epoch tine vith its corresponding 


orbital elements and the desired epoch time. This difference is 


expressed in hours. This variable will also update the elements 
to correspond to the desired epoch time. This variable normally 
should be 0.000 with the decimal point in column 7 which does 
not update the time or elements. This method of updating the 
epoch should only be used when desired epoch times and orbital 
elements are not available. 

MXNUTE®The DELTA T for time stepping to test for a moon conflict 
(in minutes). 


Card #5 

(Start date and time for this 

program) 

COL 

FORMAT 

VARIABLE 

DESCRIPTION 

1-U 

IU 

IDATE(l) 

YEAR (e.g. 1972) 

5-6 

12 

IDATE(2) 

MONTH (e.g. 12) 

7-8 

12 

IDATE(3) 

DAY (e.g. lU) 

9 

blank 



10-11 

12 

IDATE(U) 

HOUR (e.g. 02) 

12-13 

12 

IDATE(5) 

MINUTE(e.g. 35) 

1U-15 

12 

idate(6) 

SECOND (e.g. 00) 

(The seconds are 

not used) 
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Card #6 (End date and tine for this program) 


COL 

FORMAT 

VARIABLE 

description 

1-U 

IU 

IEND(l) 

YEAR (e.g. 1972) 

5-6 

12 

IEND(2) 

MONTH (e.g. 12) 

7-8 

12 

IEND(3) 

DAY (e.g. 31) 

9 

blank 



10-11 

12 

IEND(U) 

HOUR (e.g. 23) 

12-13 

12 

IEND(5) 

MINUTE (e.g. 59) 

lk-15 

12 

iend(6) 

SECOND (e.g. 00) 


(The seconds are not used) 

3.2. U JCL After Data 

If nore than one copy of the listing is desired, the Utility 
Program PATRICK may be used. The setup to inyoke PATRICK is 
as follows: 

//STEP EXEC PCM-PATRICK ,PARM* ' 9TN ,001 ,001 ' ,REOION«100K 

//INI DD DSN**.B.GO.FT09F001,DISP*(OLD,PASS) 

(This card references STEPB in the EXEC card at the beginning 
of JCL. If output is to be written to unit 10 instead of unit 
09, then FT09 in the above step must be changed to FT10.) 

NOTE: The JT09 definition must be changed if it is desired to write to 
disk rather than SYS0UT=A. The cards should be: 

//GO.FT09F001 DD DSN=i ADUM2 ,DCB= ( RECFM»VBA , LRECL*137 ,BLKSIZE*1265 ) , 

// DISP- ( NEW .PASS ) ,SPACE»( TRK , ( 50 ,10 ) ) ,UNIT=DISK 

(If FT10 is used, than change the statement to FT10F001) 

/ /0UT1 DD SYS0UT»A .SPACE* ( CYL, (5,1)), 

// DCB*(RECFM»VBA,LREC1>137,BLKSIZE«7265) 

These six cards may be repeated as many times as output is needed. 

Each set is cm exact duplicate of the other. 
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/•(End of file) 

// (Null) 

This completes the setup for the cards to run the Moon Conflicts 
Program. 
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JCL DECK FOR MOON CONFLICT PROGRAM 


//ZNSRMMON JOB ( GH650160 IB ,P ,G00392 , 004004 ), AAA ,MSGLEVEL= 1 
//B EXEC LOADER, PARM='EP=MA IN' ,REGI0N=250K 

//SYSLIN DD UNIT = (2400-9,, DEFER), LABEL= < 2 , BLP > , DI SP= < OL D, KEEP), 

/V -DCB=iRECF-M^FB-,LRECL=80 ,B LKS I ZE=3200, DEN=2 ) , VOL=SER=33462 J 
//G0.FT28F001 DD DSNAME=GDRT AP, UN I T=2400-9 , LABEL = < , BLP > , 

// VOLUME=SER=30278K,DISP=<NEW,PASS> ,DCB=(,LRECL=8276, 

// BLKSI ZE=8 28 O, D EN=2 » R ECFM = V ) 

//GO.FT09F001 DD SYSOUT=A ,D CB= ( RECFM=VBA ,LRECL=1 37, BLKSI ZE=7265 > 
//GO. FT 10F001 DD DUMMY 
• / / GO, D-ATA-5 DD- * 

130.000 130.000 1 .375 ITOS-D 

19721011 181918 

7839.840 0.00025 101.760 265.920 78.400 332.556 

0.000 50 

19721011 171800 
19721031 235000 
/* 

// 


FIGURE - 9 



3.3 Flowchart 


A flowchart of the general functions of the Moon Conflict Program is 
given in Figure 10. This diagram shows the general data flow of the 

system and the method used for obtaining satellite, sun and moon 
positions. 
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FLOW DIAGRAM OF MOON CONFLICT PROGRAM 



FIGURE - 10A 
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FLOW DIAGRAM OF MOON CONFLICT PROGRAM (CONT’D) 



FIGURE - Iob 
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3 . 4 Subroutines 

Below are listed all the subroutines necessary in the Moon Conflict 
Program. All of these routines are mentioned or described in the 
Subroutine section of the DSAS Illumination Program description. 
Refer to that section (2.4) for further details. 


BRHV 

DMAD 

BBHWR 

ELRV 

BRWR1 

MAD 

BRWR2 

XKEP 

ALLOT 

RQ4M2T 

ALLOTZ 

RDJPL 

ATAHQ 

DELTAD 

DJUL 


Output Description 

(see Figure 11 ) 


A printed listing is the only output from the Moon Conflict Program. If 
tape output is desired, the JCL must be changed to write a tape. A 
description of the output is as follows: 

1. The first page is a title page containing the program i dentifi cation - 

"Moon Conflict Test for * Satellite. 

* ftie satellite name will be entered in this space. 

2. The second page contains pertinent information about the mounting 
angles of the horizon scanners from the positive spin axis, the 
epoch time in year, month, day, hours, minutes and seconds, and 

the orbital elements used in this program* Also included is the start 
and end times of the program. 
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3. Hie next page and any following pagea contain the times during a 
Booth vhen each scanner will be in conflict with the noon. Each 
time is given in year, month , day, hour, and minute , and corresponding 
to each time is given the right ascension in degrees of the sun and 
moon. The column labelled "difference" gives the difference in degrees 
betveen the right ascension of the sun and noon on the basis of 180 
degrees. When the difference is greater than 90 degrees, then an 
actual conflict will occur. Vhen the difference is less than 90 degrees, 
then the illuminated disk of the moon is not enough to cause a real 
conflict with the horizon scanner. 
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